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With the aim of constraining the structural properties of 
neutron stars and the equation of state of dense matter, we 
study sudden spin-ups, glitches, occurring in the Vela pulsar 
and in six other pulsars. We present evidence that glitches 
represent a self-regulating instability for which the star pre- 
pares over a waiting time. The angular momentum require- 
ments of glitches in Vela indicate that > 1.4% of the star's 
moment of inertia drives these events. If glitches originate 
in the liquid of the inner crust, Vela's 'radiation radius' Rac 
must exceed ~ 12 km for a mass of I.AMq. Observational 
tests of whether other neutron stars obey this constraint will 
be possible in the near future. 



I. INTRODUCTION 

The sudden spin jumps, or glitches, commonly seen in 
isolated neutron stars are thought to represent angular 
momentum transfer between the crust and the liquid in- 
terior In this picture, as a neutron star's crust spins 
down under magnetic torque, differential rotation devel- 
ops between the stellar crust and a portion of the liquid 
interior. The more rapidly rotating component then acts 
as an angular momentum reservoir which occasionally ex- 
erts a spin-up torque on the crust as a consequence of an 
instability. The Vela pulsar, one of the most active glitch- 
ing pulsars, typically undergoes fractional changes in ro- 
tation rate of ^ 10~^ every three years on average 
With the Vela pulsar having exhibited 13 glitches, mean- 
ingful study of the statistical properties of these events 
is now possible. 

In this Letter we study the time distribution of Vela's 
glitches and determine the average angular momentum 
transfer rate in Vela and in six other pulsars. We present 
evidence that glitches in Vela represent a self-regulating 
instability for which the star prepares over a waiting in- 
terval. We obtain a lower limit on the fraction of the 
star's liquid interior responsible for glitches. Assuming 
that glitches are driven by the liquid residing in the inner 
crust, as in most glitch models, we show that Vela's 'radi- 
ation radius' is i?oo ^12 km for a mass of \.AMq. Future 



measurements of neutron star radii will check the univer- 
sality of this constraint and hence test our understanding 
of neutron star structure and the origin of glitches. 



II. REGULARITY OF ANGULAR MOMENTUM 
TRANSFER 

A glitch of magnitude Af7i requires angular momen- 
tum 



A Ji = IcAVti 



(1) 



where 1^ is the moment of inertia of the solid crust plus 
any portions of the star tightly coupled to it. Most of 
the core liquid is expected to couple tightly to the star's 
solid component, so that Ic makes up at least 90% of the 
star's total moment of inertia ||^. Glitches are driven by 
the portion of the liquid interior that is differentially ro- 
tating with respect to the crust. The cumulative angular 
momentum imparted to the crust over time is 
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where f2 — 70.4 rad is the average spin rate of the 
crust over the period of observations. Fig. 1 shows the 
cumulative dimensionless angular momentum, J{t)/Icn, 
over 30 years of glitch observations of the Vela pulsar, 
with a linear least-squares fit. The average rate of angu- 
lar momentum transfer associated with glitches is Ic^A, 
where A is the slope of the straight line in Fig. 1: 



(6.44 ±0.19) X 10"' yr 



(3) 



This rate A is often referred to as the pulsar activity 
parameter. 

The angular momentum flow is extremely regular; 
none of Vela's 13 glitches caused the cumulative angu- 
lar momentum curve to deviate from the linear fit shown 
in Fig. 1 by more than 12%. To assess the likelihood 
that the linear trend could have arisen by chance, we 
tested the statistical robustness of this result. We gener- 
ated many sets of simulated data in which the occurrence 
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times of the glitches remained as observed, but the mag- 
nitudes of the 13 glitches were randomly shuffled. We 
compared the observed to those for the deviations of 
the randomly shuffled data from linear fits. The for 
the shuffled data was less than that of the real in 
only ^ 1.4% of cases, strongly suggesting that the rate of 
angular momentum flow associated with glitches is rea- 
sonably constant. 

Additionally, the near uniformity of the intervals be- 
tween the glitches in Fig. 1 suggests that glitches occur 
at fairly regular time intervals. The standard deviation 
in observed glitch intervals is 0.53(At), where (At) = 840 
d is the average glitch time interval. The probability of 
13 randomly-spaced (Poisson) events having less than the 
observed standard deviation is only ~ 1%. 

The data of Fig. 1 indicate that Vela's glitches are 
not random, but represent a self-regulating process which 
gives a relatively constant flow of angular momentum to 
the crust with glitches occurring at fairly regular time 
intervals. 



III. THE GLITCH RESERVOIR'S MOMENT OF 
INERTIA 

The average rate of angular momentum transfer in 
Vela's glitches constrains the properties of the angular 
momentum reservoir that drives the spin jumps. In par- 
ticular, the frequent occurrence of large glitches requires 
that a significant fraction of the interior superfluid spins 
at a higher rate than the crust of the star. Between 
glitches, the reservoir acquires excess angular momentum 
as the rest of the star slows under the magnetic braking 
torque acting on the crust. Excess angular momentum 
accumulates at the maximum rate if the reservoir does 
not spin down between glitches. Hence, the rate at which 
the reservoir accumulates angular momentum capable of 
driving glitches is limited by 



A similar analysis for six other pulsars yields the results 
shown in Fig. 2. An earlier analysis of glitches in Vela 
gave Ifcs/Ic > 0.8% After Vela, the most significant 
limit is obtained from PSR 1737-30 which gives Ircs/Ic > 

Gl737 = !%■ 

The similarity of G for the five objects of intermediate 
age suggests that glitches in all these objects are driven 
by internal components with about the same fractional 
moment of inertia. In terms of G, the Crab pulsar and 
PSR 0525-1-21 appear to be unusual. It may be that 
the Crab's angular momentum reservoir loses its excess 
angular momentum between glitches, perhaps through 
thermal creep of superfluid vortices (see, e.g., [^). The 
value of G for PSR 0525-1-21 is not well determined, since 
only two glitches from this object have been measured. 



IV. IMPLICATIONS FOR THE DENSE MATTER 
EQUATION OF STATE 

The constraint of I^es/Ic > 1.4% for Vela applies re- 
gardless of where in the star glitches originate. Many 
glitch models, however, assume that the internal angular 
momentum reservoir is the superfluid that coexists with 
the inner crust lattice where the pinning of super- 
fluid vortex lines sustains a velocity difference between 
the superfluid and the crust. Here we explore the im- 
plications of this interpretation. We begin by describing 
how the moment of inertia of the superfluid in the neu- 
tron star crust relates to the nuclear matter equation 
of state (EOS) and the observable properties of neutron 
stars. 

Ravenhall & Pethick have shown that, for various 
equations of state, the total moment of inertia / is given 
by the approximate expression 
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J res ^ -^rcsl^l 



(4) 



where (l is the average spin-down rate of the crust, and 
/res is the moment of inertia of the angular momentum 
reservoir (not necessarily one region of the star). Equat- 
ing Jios to the average rate of angular momentum transfer 
to the crust, iS^A, gives the constraint. 



b^>^A = G, 

Ic ~ 



(5) 



where the coupling parameter G is the minimum fraction 
of the star's moment of inertia that stores angular mo- 
mentum and imparts it to the crust in glitches. Using 
the observed value of Vela's activity parameter A and 
f2/|Ji| = 22.6 Kyr, we obtain the constraint 



where p is the mass-energy density, P is the pressure, 
and e"^ is the local gravitational redshift. This expres- 
sion, which holds in the limit of slow rotation, deflnes 
the integral J. This integral can be evaluated following 
Lattimer & Prakash |^ who noted that p (x 1 — (r/i?)^ 
throughout most of the interior of a neutron star (but not 
in the crust), for all commonly- used equations of state. 
With this approximation, it can be shown |^ that 



J~ -MR^A, 
7 

.2\-l 



(8) 



where A = (1 - 2GM/Rc^)-^ and M is the total steflar 
mass. 

Equation can also be used to determine the moment 
of inertia of the crust plus liquid component. In the crust 
P is << pc^, and the TOV equation is 
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Using this approximation in eq. (|^) gives the fraction of 
the star's moment of inertia contained in the sohd crust 
(and the neutron hquid that coexists with it): 
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Here Ai? is the radial extent of the crust and Pt is the 
pressure at the crust-core interface. A similar approxi- 
mation is obtained in Ref. 7 (equation 17); either ap- 
proximation is adequate for the estimates we are making 
here. In most of the crust, the equation of state has the 
approximately polytropic form P oc p^/"^ , giving H 
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8Pt 4.5+(A-l)~ 
ntmnC^ A — 1 



(11) 



where rit is the density at the core-crust transition and 
nin is the neutron mass. A/// can thus be expressed as a 
function of M and R with an additional dependence upon 
the EOS arising through the values of Pt and rit- How- 
ever, Pt is the main EOS parameter as nt enters chiefly 
via a correction term. In general, the EOS parameter Pt 
varies over the range 0.25 < Pt < 0.65 MeV fm~'^ for 
realistic equations of state |^ . Larger values of Pt give 
larger values for A///, as can be seen from eq. [|ll|. 

Combining of eqs. ||T^ and with a lower limit on 
AI and an upper limit on Pt gives a lower limit on the 
neutron star radius for a given mass. In order to relate 
our observational bound on I-^cs/Ic to AI, we assume 
that the angular momentum reservoir is confined to the 
neutron superfluid that coexists with the nuclei of the 
inner crust. In this case, I^cs < AI and Ic > I — AI. 
Our observational limit on I^^s then gives AI / {I — AI) > 
AI/Ic > Iics/Ic > 0.014. To obtain a strong lower limit 
on the neutron star radius, we take Pt — 0.65 MeV fm^"^ 
and nt = 0.075 fm~'^. Combining the relations |Jo| and 
pl| , gives the heavy dashed curve in Fig. 3. This curve 
is given approximately by 



i? = 3.6 -f 3.9M/M0 . 



(12) 



Stellar models that are compatible with the lower bound 
on /res must fall below this line. Smaller Pt reduces the 
crustal moment of inertia and gives a more restrictive 
constraint. For example, Pt = 0.25 MeV fm^"^ moves 
the constraining contour to approximately R = 4.7 -I- 
4.IM/M0 (thin dashed curve of Fig. 3). 



V. DISCUSSION 

To summarize our conclusions regarding the statistics 
of Vela's glitches, we find that angular momentum is im- 
parted to the crust at regular time intervals at a rate that 
has remained nearly constant for ~ 30 yr. These data 
narrowly constrain the coupling parameter G which is the 



minimum fraction of the star's moment of inertia that is 
responsible for glitches. For Vela we find G = 0.014, in- 
dicating that least 1.4% of the star's moment of inertia 
acts as an angular momentum reservoir for driving the 
glitches, regardless of where in the star this angular mo- 
mentum reservoir is, or how it is coupled to the crust. 
Variation of G by a factor of less than ^ 3 for stars in 
the age group 10* — 10^ yr suggests that glitches in stars 
in this age group all involve regions of about the same 
fractional moment of inertia. 

Mass measurements of radio pulsars in binary sys- 
tems and of neutron star companions of radio pulsars 
give neutron star masses consistent with a very nar- 
row distribution, M = 1.35 ± O.O4M0 indicated by 
the pair of horizontal dotted lines in Fig. 3. If Vela's 
mass falls in this range, eq. |^ constrains i?^8.9 km, 
under the assumption that glitches arise in the inner 
crust superfluid. The quantity constrained by observa- 
tions of the stellar luminosity and spectrum is the 'ra- 
diation radius' R^ = A^^^R = (1 - 2GM/Rc^)'^/^R. 
If M = 1.35Af0 for Vela, the above constraint gives 
i?oo <^ 12 km if glitches arise in the inner crust. For com- 
parison, we show in Fig. 3 the mass-radius curves for sev- 
eral representative equations of state (heavy solid lines). 
Measurement of Roo ^ 13 km would be inconsistent with 
most equations of state if M ~ 1.35M0. Stronger con- 
straints could be obtained if improved calculations of nu- 
clear matter properties indicate Pt significantly less than 
0.65 MeV fm-3. For example, for M ~ 1.35M0, i?oo ^ 13 
km would be required if Pt = 0.25 MeV fm""^. A mea- 
surement of Roo ^11 km would rule out most equations 
of state regardless of mass or the angular momentum re- 
quirements of glitches. 

A promising candidate for a decisive measurement of 
a neutron star's radiation radius is RX J185635-3754, an 
isolated, non-pulsing neutron star ||l^. A black body 
fit to the X-ray spectrum gives i?oo = 7.3(D/120 pc) 
km where D is the distance (known to be less than 120 
pc). However, either a non- uniform surface tempera- 
ture or radiative transfer effects in the stellar atmosphere 
could raise this estimate significantly |^^. HST obser- 
vations planned for this year should determine the star's 
proper motion and parallax, and hence, the distance. Fu- 
ture CHANDRA observations should yield more detailed 
spectral data and could establish the composition of the 
atmosphere if absorption lines are identified. These dis- 
tance and spectral data may establish whether this ob- 
ject's radius is consistent with an inner crust explanation 
of neutron star glitches. 
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Fig. 1 - Cumulative dimensionless angular momentum, 

J/Ic^, imparted to the Vela pulsar's crust as a function 
of time. The straight line is a least-squares fit. 
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Fig. 2 The coupling parameter G. The strongest con- 
straints are obtained for Vela and PSR 1737-30, for which 
13 and 9 glitches have been observed, respectively. Dia- 
monds indicate objects with only two observed glitches, 
for which error bars could not be obtained. References; 
0525H-21 [14], Crab [15], Vela [2], 1338-62 [16], 1737-30 
[17,18], 1823-13 [18]. 
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Fig. 3 - Limits on the Vela pulsar's radius. The heavy 

dashed curve delimits allowed masses and radii that are 
compatible with the glitch constraint A//(/ — A/) > 
1.4% for Pt = 0.65 MeV fm^'l The thin dashed curve 
corresponds to Pt = 0.25 MeV fm"-^ and gives a more 
stringent though less conservative constraint. The dot- 
dashed curve corresponds to A//(J — A/) > 2.8% and 
Pt = 0.65 MeV. The horizontal dashed lines indicate the 
mass limits for the survey of 26 radio pulsars of Ref. 19. 
Also displayed are mass-radius relations for the equations 
of state of Akmal & Pandharipande [10] (curves a and b), 
Wiringa, Fiks & Fabrocini [11] (curves c and d), Miiller & 
Serot [12] (curves e and f) and the kaon EOS of Glenden- 
ning & Schaffner-Bielich [9] (curves g and h) . The crosses 
indicate where a given EOS has A//(/ - A7) = 1.4% 
(with Pt = 0.65 MeV fm~^). Curves without crosses 
have AI/{I - A/) > 1.4% for all stable R. Thin curves 
are contours of constant radiation radius Roo ■ 
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